Magnetic phase diagrams of barcode-type nanostructures 
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The magnetic configurations of barcode-type magnetic nanostructures consisting of alternate fer- 
romagnetic and nonmagnetic layers arranged within a multilayer nanotube structure are investigated 
as a function of their geometry. Based on a continuum approach we have obtained analytical ex- 
pressions for the energy which lead us to obtain phase diagrams giving the relative stability of 
characteristic internal magnetic configurations of the barcode-type nanostructures. 
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I. INTRODUCTION 



Magnetic nanoparticles have attracted increasing in- 
terest among researchers of various fields due to their 
promising applications in hard disk drives, magnetic ran- 
dom access memory, and other spintronic devices. [H-IH 
In addition, these magnetic nanoparticles can be used for 
potential biomedical applications, such as magnetic res- 
onance imaging (the nanoparticles can be used to trace 
bioanalytes in the body), cell and DNA separation, and 
drug delivery. Q To apply nanoparticles in various po- 
tential devices and architectures, it is very important to 
control their size and shape in order to keep the thermal 
and chemical stability. [7| 

The trusty sphere remains the preferred shape for 
nanoparticles but this geometry leaves only one surface 
for modification, complicating the generation of multi- 
functional particles. Thus, a technology that could mod- 
ify differentially the inner and outer surfaces would be 
highly desirable. Q Tubular nanostructures have stim- 
ulated extensive research efforts in recent years because 
of their particular significance for prospective applica- 
tions. A wide range of materials including semicon- 
ductors, polymers, and metals have been prepared in 
the form of nanotubes. [9l4l3l| Although the magnetic 
nanotubes has been intensely investigated, barcode-type 
nanostructures have received less attention, in spite of 
tailoring the multisegmented nanotube structure, along 
with the functionalization of the inner wall surface of 
barcode-type nanotubes with various molecules (for ex- 
ample, proteins and DNA). Moreover, they are expected 
to be particularly useful in the field of catalysis, ad- 
vanced microfiuidics, molecule separation and biological 
and magnetic sensors as well. It is worth to 

mention that barcode-type magnetic nanostructures con- 
sisting of regular arrays of magnetic segments have been 
considered as providing the basis for extending magnetic 
storage densities beyond the superparamagnetic limit. 
In such system, a single tube with n magnetic layers 
might store up to 2 n bits, whose volume is much larger 
than those of the grains in conventional recording media, 



bearing this way thermal fluctuations and increasing the 
recording density by a factor 2 n_1 . [13, [HI Recently, [l4| 
the preparation of metallic nanotubes based on the pref- 
erential electrodeposition of a metal along the pore walls 
of an anodic alumina oxide (AAO) membrane, in the 
presence of metallic nanoparticles on the wall surfaces, 
has been reported. In the paper by Lee et al. [l4j they 
were able to prepare multisegmented metallic nanotubes 
with a bimetallic stacking configuration along the tube 
axis, showing different magnetic behavior as compared 
with continuous ones, which encourage a study about 
the possible magnetic configurations and magnetostatic 
interactions in these barcode-type magnetic nanotubes. 
Clearly, for the development of magnetic devices based on 
those arrays, knowledge of the internal magnetic struc- 
ture of the barcode-type nanostructures is of fundamental 
importance. 

The purpose of this paper is to investigate the magnetic 
ordering of barcode-type nanostructures. Our particles 
are characterized by a set of geometrical parameters, as 
depicted in Fig. 1. First, their external and internal 
radii, R and a, respectively, and total length, L, which 
includes the magnetic matter as well as the nonmagnetic 
one. It is convenient to define the ratio ft = a/R, so that 
ft — represents a solid cylinder and ft — > 1 corresponds 
to a very narrow tube. We denote with W the length of 
each ferromagnetic segment and with d the length of each 
nonmagnetic portion, so that, if we have a barcode-type 
structure with n magnetic tubes, the total length can be 
written as L — nW + (n + 1) d. 



II. MODEL AND DISCUSSION 

We adopt a simplified description of the magnetic sys- 
tem in which the discrete distribution of magnetic mo- 
ments is replaced by a continuous one, defined by the 
magnetization vector field M(r) such that M(r)5w gives 
the total magnetic moment within an elementary vol- 
ume Sv centered at r. The total magnetic energy (E to t) 
is generally given by the sum of four terms: exchange, 
dipolar, anisotropy and Zeeman contributions, which are 
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FIG. 1: Geometrical parameters and magnetic configurations 
of barcode-type nanostructures. 



taken from the well known continuum theory of ferro- 
magnetism. [22[ As we are interested in the study of 
the relative stability of the zero-field magnetic ground 
states, the contribution of the Zeeman energy can be dis- 
regarded. Under these assumptions, the magnetic energy 
is just given by the dipolar (E dip ), exchange (E ex ), and 
anisotropy (Ek) contributions. 

The total magnetization can be written as M (r) = 
Ei=i Mj (r), where Mi (r) is the magnetization of the 
i-th ferromagnetic segment. In this case, the magneto- 
static potential U (r) splits up into n components, f/j (r), 
associated with the magnetization of each ferromagnetic 
segment. Then, the total dipolar energy can be written 
as E dip = Yh=i E di P (*)+Z)r=i Y,j=i+i E int where 

Edip (*) = y J M * ( r ) ' Vf/ * M du ' 

is the dipolar contribution to the self-energy of the i-th. 
ferromagnetic segment, and 



E, 



u (*, j) = A*o / M * (') • VET,- (r) dv, 



is the dipolar interaction between ferromagnetic seg- 
ments i and j. 

Usually, the exchange energy E ex in multilayer nanos- 
tructures has contributions from both, the direct ex- 
change interaction within the magnetic segments and the 
other from the indirect interaction between them medi- 
ated by the conduction electrons in the nonmagnetic lay- 
ers. Since the indirect interaction decays rapidly with the 
thickness of the nonmagnetic segment, it can be neglected 
provided d is large enough. A good estimate of the range 
of the indirect exchange interaction can be obtained from 
the results for multilayers. [23| As a general result, the in- 
terlayer exchange coupling vanishes for spacer thicknesses 
greater than a few nanometers, which does not exceed the 



value of the exchange length l x = y/2A/ hqMq of ferro- 
magnetic metals. Here we focus our attention on those 
cases in which d is not smaller than the magnetic mate- 
rial's l x , as the tubes fabricated by Lee et al. [IH which 
satisfy d 3> l x and thus interlayer exchange coupling can 
be safety neglected. Therefore, to a good approximation 
can be written as E ex — Y^—y E ex (i), where E ex (i) = 

(Vrrii Z ) 2 dv. Here, rrii = 



A J (Vra«) + (Vm 



HI J 



s ) = Mi /Mo is the magnetization normal- 
ized to the saturation magnetization Mq and A is the 
stiffness constant of the magnetic material. 

The cubic anisotropy energy of the particle can be 
added by means of the following expression: 

and the uniaxial anisotropy energy is given by 

E u (i) = ~K U J mf z dv. 

On the basis of the above results, the total energy 
of the barcode-type nanostructure can be written as 
Etot = J2i=i E self(i) + Y%=lT%=i + i E int(i,j), where 
E se i f (i) = E dip (i) + E ex (i) + E k (i) is the self-energy of 
the ferromagnetic segment i, and Ei nt is the (dipolar) 
interaction energy between two magnetic segments. We 
will proceed to describe the magnetization of the different 
states we are considering here and then we will evaluate 
the magnetic energy of each configuration. Results will 
be given in units of fioM§l x , i.e., E = E/^ MqI x . 

A. Magnetic configurations 

It has been shown recently that single magnetic nanor- 
ings present three basic ground states depending on their 
geometry (see Fig. 1). (13, H^j These configurations are: 
(Fz) a quasi uniform magnetization state oriented in the 
direction parallel to the cylindrical axis (z axis); (Fx) a 
quasi uniform magnetization state oriented in the plane 
perpendicular to the z axis; and (V) a flux-closure vor- 
tex state. For long nanorings (W ^> R), the Fx phase is 
not present, [24l . [25| a result that holds for nanotubes. 
[2ll22l 

It has been shown for Rothman et al. [28| that for 
magnetized nanorings the single-domain in-plane ground 
state is the onion state. Besides, the single-domain axial 
state for rings with small inner diameter might be similar 
to the flower state expected in thick axially magnetized 
cylinders. [2^] Therefore, it may appear questionable to 
select the uniform in-plane and axial as single-domain 
states to build the phase diagram upon, as they are not 
stable configurations at a zero applied field. However, it 
has been verified by micromagnetic simulations [2~ij and 
analytical calculations [25| that the energy difference be- 
tween the actual single-domain ground state in a nanor- 
ing, and the uniform state, often turns out to be very 
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small. From this available evidence, we conclude that 
the replacement of the more correct quasi-uniform states 
by simpler ideal uniform states, only results in uncertain- 
ties on the exact location of the phase boundaries and on 
some physical values extracted from the phase diagram. 



1. Fz state 

For the Fz state, where the magnetization of the n 
ferromagnetic segments is uniform and parallel to the 
z axis, M (r) can be approximated by MqZ, where z is 
the unit vector parallel to the axis of the nanotube. In 
this case the exchange contribution to the self energy 
vanishes, and the reduced self energy takes the form [261 ] 
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FIG. 2: Stray field direction (arrows) generated by a single 
nanotube (left picture) and a multisegmented nanotube (right 
picture) magnetized in the +z direction. These two examples 
correspond to nanotubes with the same geometrical parame- 
ters as those studied experimentally by Lee et al 



where J\ (z) is a Bessel function of the first kind and 
k u = 2K U / hoMq . In order to calculate the interaction 
energy between the ferromagnetic segments, we first need 
to calculate the magnetostatic potential U (r) of a single 
tubular structure. The expression for this potential has 
been previously reported [30[ and is given by 



71 ,T pOO 11 

U{r,z) = -£ ^ Jo (At) [RJ x (kR) - aJ x (ka)] 



-+z 



From this equation it is possible to obtain the expression 
for the magnetostatic field. Thus we write, H (r, z) = 
-S7U (r, z) = H r (r, z)f + H z (r, z) z with 

H r (r, z) = — - j dk Ji (kr) [RJ X (kR) - aJ x (ka)] 



al. [2l|, l30( and is given by 
ttR 3 f°° dq 



El\ [d] 



?3 / a- 
l x Jo Q 



-e H R 



(l - e-«$) (Ji (q) - 0Ji (qf3)f 



Thus, the reduced total energy for the Fz state can be 
expressed as 



E tot = nE sIlf 



where g z (n, q, a) 



ttR- 



f°° dq „w / 

Jo r v 



( Ji (q) - PJi (qfi)) 9z(n,q,<r), 



fn-l)e*^+e 



2. Fx state 



and a 



d+W 
R ■ 



and 



H z (r, z) 



where Y (W, z) 



For the Fx state, M (r) can be generally considered as 
M cos [(i — 1) 9} x+M sin [(£ — 1) 9] y, which represent a 
helicoidal magnetic state, with 9 the angle between the 
in-plane magnetization of adjacent segments. For the in- 
dk J (kr) [RJ X (kR) - aJ x (ka)} Y (W, z) plane state, the exchange and anisotropy contributions to 

the self energy vanish and the reduced self energy takes 
the form [Hf 



sign {W_ z } e -k\^-*\ _ 

sign (— — z) e~ k \~~~~ z \ . The function sign(x) 
gives —1, or 1 depending on whether x is negative, 
zero, or positive. Figure 2 illustrates the magnetostatic 
field profile calculated analytically for nanotubes with 
the same geometrical parameters as the ones investigated 
experimentally by Lee et al. [lij]. Finally, the reduced 
interaction energy between two tubular nanostructures 
has been calculated in rather general way by Escrig et 
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(Jl (q) - PJi (q(i)Y 



where J\ (z) is a Bessel function of the first kind. In order 
to calculate the interaction energy between the ferromag- 
netic segments, we first need to calculate the magneto- 
static potential U (r) of a single tubular structure. The 
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expression for this potential is given by 

dk 



U (r,<f),z) = M cos <j) 



o k 



Ji{kr)f{k) 



e kz sinh (k^- 
i~ k -r cosh (kz 
e kz sinh (k f ) 



z > — 

-f <z<f 



z < 
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where / (k) = (RJ\ (kR) — aJ\ (kaj). Finally, the re- 
duced interaction energy between two tubular nanostruc- 
tures is given by 



TrR 3 

~w 



dq 



l-e-«TT (Ji (q) - fiJi (q/3)Y 



Thus, the reduced total energy for the Fx state can be 
expressed as 



E fot = nE self 



(Ji(q)-PJi(qP)fg* (n,q,a,e), 
where g x (n, q, a, 9) = ELi+i e-^-i) cos [(j - i) 0] 



and a — . From this expression we can conclude 
that for zero applied field the total energy of this state 
is further reduced for the value 9 — ir, independently of 
the value of n. Thus, for the particular case when 9 = tt 



we obtain g x (n, q, a, tt) 
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FIG. 3: Phase diagrams for barcode-type nanostructures giv- 
ing the regions in the RW plane where one of the configura- 
tions has lower energy. We have used P = 0.5, d = l x , and 
n = 5. 



TABLE I: Parameters for different materials taken from Ref. 
[31]. Uniaxial (K u ) cobalt is denoted with a superscript *. 
Iron, permalloy and nickel have cubic (K c ) anisotropy. 



Material 


K (J/m 3 ) 


K 


Iron 


4.8 x 10 4 


0.0264 


Cobalt* 


4.1 x 10 5 


0.3329 


Permalloy 


-3.0 x 10 2 


-0.0007 


Nickel 


-4.5 x 10 3 


-0.0304 



3. Vortex state 

Finally, for the vortex state V, M (r) can be approx- 
imated by Mo0, where is the azimuthal unit vector. 
Due to the condition of perfect flux closure in the vortex 
configuration, one magnetic nanostructure in such con- 
figuration does not interact with others, independently 
of the magnetic configuration of those. Thus, there is no 
difference between clockwise and counter-clockwise di- 
rections. Finally, the reduced total energy for the vortex 
state is given just by the n self energies |2l|, [26| 



h = — n 



k c nWR 2 . 



P 2 



Here, k c = 2K c /fi M§. 



B. Phase diagram for multisegmented nanorings 

We proceed to investigate the relative stability of the 
configurations. Phase diagrams are shown in Fig. 3 for 
d — l x , P = 0.5, and n = 5. Anisotropy for four different 
materials are considered according to values presented in 



Table 1. The diagrams show three regions, corresponding 
to configurations Fz, Fx, and V, as in the case of a single 
nanoring (n = I). Notice that for the case of Co, the 
existence of a strong uniaxial anisotropy favors the Fz 
phase, decreasing the other two phases, specially the V 
one. In the case of a cubic anisotropy, the transition 
lines are similar to the case of a phase diagram without 
anisotropy. Because of its very low anisotropy, results 
for permalloy describe reasonably well a material with 
no anisotropy, as it was pointed in Ref. [27]. 

Since nanostructures are usually polycrystalline, the 
crystallographic orientations of the crystallites are ran- 
dom and, as a consequence, the average magnetic 
anisotropy of the particle is very small. In view of that, 
it will be neglected in our calculations. [HI, [33| 

For different values of n we can determine the ranges of 
values of the dimensionless radius R/l x and length W/l x 
within which one of the three configurations is of lowest 
energy. The boundary line between any two configura- 
tions can be obtained by equating the expressions for the 
corresponding total energies. Figure 4 illustrates phase 
diagrams for d = l x , (3 = 0.5, and n = 1 (solid lines), 3 
(dotted lines) , and 5 (dashed lines) . It is important to ob- 
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serve that for the Fz and Fx states the exchange energy is 
the same. Then, in the absence of applied magnetic fields 
and crystalline anisotropies, the dipolar energy is funda- 
mental to obtain the magnetic configuration of lowest en- 
ergy. Thus, the dipolar contribution represents the shape 
anisotropy that, for multisegmented nanostructures with 
a small length (namely nanorings) the low energy state is 
the quasi uniform in-plane configuration Fx. [25| As the 
length is increased, but keeping the radius small enough, 
there is a transition to the out-of-plane state Fz at a 
critical length whose value depends on R, /?, and the 
exchange length l x . As the radius is increased, the mag- 
netizations turns to the vortex configuration at a critical 
radius depending on the values of W, (3, and l x . Finally, 
by comparing our results we observe differences in the 
behavior of the triple point as a function of n. The triple 
point occurs for smaller R/l x when n is decreased. 




FIG. 4: Phase diagrams for barcode-type nanostructures with 
P — 0.5 and d = l x . 

Similar to the case of a single ring, the phase diagram 
changes with j3 (25|. The dependence of the whole dia- 
gram on the value of n can be investigated by looking 
at the trajectories of the triple point in the RW plane 
as functions of (3. Such trajectories are shown in Fig. 
5 for d = l x and different values of n. We remark that 
the radius Rt of the triple point represents the smallest 
value of R for which the vortex configurations are stable, 
and Wt is the biggest value of W for which the in-plane 
configurations are stable. 



C. Phase diagram for multisegmented nanotubes 

As the multisegmented tubes that motivate this work 
[Tij satisfy W/R S> 1, then the Fx phase can be left 
out of consideration. Thus, to obtain an expression for 
the transition line separating the Fz phase from the V 
phase we match the expressions for the energy of these 
two configurations. It is important to mention that, for 
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FIG. 5: Trajectories of the triple point in the phase diagrams 
in Fig. 4 as functions of P, for n = 1 (squares), 3 (circles), 
and 5 (triangles). 



tubes with long radius, it has been observed a third state 
which is a mixture of the other two and has been called 
bamboo or mixed state. [3J-|37[ As it is known, the con- 
sideration of no uniform magnetic configurations compli- 
cates considerably the calculations and for simplicity, we 
studied multisegment magnetic nanotubes whose radius 
are not big enough to allows the formation of relevant 
vortex domains at the extremes of the tube. Figure 6 
presents the transition line for n = 1 and n — 2. To the 
left of each line Fz state prevails while to the right of 
the same line the vortex V configuration is more stable. 
Labelled dots (continuous) and (multisegmented) in Fig. 
6 correspond to the cases of the two hysteresis curves 
reported in the experimental paper by Lee et al. fl4| de- 
fined by (continuous) n = 1, R = 150 nm, W = 16 /im, 
/3 = 0.75, and l x — 8.225 nm; (multisegmented) n = 2, 
R = 150 nm, W = 800 nm, d = 4800 nm, (3 = 0.75, and 
l x = 8.225 nm. It is important to note that the transi- 
tion line for n = 1 is almost equal to the one with n = 2. 
It due to that average distance between the neighboring 
Ni segments was big enough (d = 4.8 ^m) avoiding thus 
the interaction between the segments. From this figure 
we can conclude that the multisegmented system is well 
inside the V phase while the continuous system is inside 
the Fz phase. It allows us to understand why the ex- 
perimental samples show a different magnetic behavior; 
simply they have substantial differences in their length 
of the ferromagnetic segments. 

The results presented above may be generalized. We 
now proceed to investigate the transition line separating 
the Fz phase from the V phase. To obtain an expres- 
sion for this transition line we match the expressions for 
the energy of these two configurations. This leads to 
W/l x = a(P) x R 3 /l x . Function a(f3) is plotted in Fig. 
7. Care must be applied in the limits of the intervals 
for j3. In particular, when ft goes to 1 we deal with ex- 
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FIG. 6: Magnetic phase diagrams of non-interacting multiseg- 
mented nanotubes for different values of n. The dimensions 
of the tube, W and R, are normalized to the exchange length 
l x . Experimental points are discussed in the text. 



tremely narrow nanotubes, where eventual surface rough- 
ness and thickness irregularities of the nanotubes become 
important. On the other side, when /3 goes to zero we 
are approaching the limit of a solid cylinder, where the 
core in the vortex phase becomes important and must be 
considered to get the solution. As the multisegmented 
nanotubes considered experimentally have (3 ks 0.75, we 
have neglected these two cases. 



lines for the production of nanostructures with techno- 
logical purpose. 
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III. CONCLUSIONS 

In conclusion, we have studied the relative stability 
of ideal configurations of magnetic barcode-type tubular 
nanostructures composed of alternate ferromagnetic and 
non-magnetic layers. In such systems we investigated the 
size range of the geometric parameters for which different 
configurations are of lowest energy. Results are summa- 
rized in phase diagrams which clearly indicate that the 
magnetic behavior of such structures can be tailored to 
meet specific requirements provide a judicious choice of 
such parameters is made. The lines separating the mag- 
netic phases and, in particular, the triple point, are very 
sensitive to the geometry of the barcode-type nanostruc- 
tures. The phase diagrams presented can provide guide- 
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FIG. 7: Function ot{0) defining the transition condition from 
the phase diagram of a multisegmented nanotube. 
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